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Abstract. Transport through a point contact is accurately modelled by assigning 
to the junction an ensemble of independent transport channels with possibly different 
transmissions. We here argue that for a series of two contacts, coherently coupled 
across an island, the transport channels are different from the ensembles that would 
describe each contact taken as stand-alone device. We further show that instead of two 
sets of channels with manifold cross-links over the island the double junction can be 
described by pairs of channels from both sides coherently coupled together, where each 
pair, however, has no coherent connection to the others. This finding will substantially 
simplify modelling transport by a Green's functions technique. Additional channels 
through only one junction may complete the picture. Finally we discuss how partial 
coherence across the island with an appropriate ansatz can be modelled in the same 
scheme. 



PACS numbers: 05.60.Gg, 02.10.Yn, 73.23.-b, 73.63.Rt 
1. Introduction 

In low-dimensional electron systems [H Ej [3] the quantization of conductance [1] is 
observed. For constrictions in two-dimensional electron gases the conductance adopts 
multiples of the conductance quantum Go = /h because the transverse wave numbers 
determine how many modes contribute [5l[6]. Single-atomic size contacts exhibit typical, 
reproducible and material dependent conductances, which, however, in general are no 
integer multiples of Go [Z]- The conductance is associated to an ensemble of transport 
channels, the number of which corresponds to the number of valence orbitals of the 
element used. The individual channel transmissions in the range ... 1 reflect the wave- 
function overlap from the central to neighbouring and following atoms |8j. Although 
atomistic ab initio calculations greatly complement experiments in this field, they are not 
needed to deduce transport channels and their transmissions from measurements. The 
description using an ensemble of transport channels is a more general concept applicable 
to any sort of point contact [9j. The contact - including leads to bulk electrodes - can be 
viewed as a black box that behaves like the associated ensemble of transport channels. 
A deeper interpretation of the channels is not required. We shortly review here, in our 
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own way, the theoretical construction behind the description of a single contact in terms 
of a channel ensemble, because this basis is needed for generalizing the concept to the 
double junction. 

2. Single junction 

Regarding the left and right side of a point contact as in Fig. la at first as uncoupled 
reservoirs let there be a finite number of orthonormal modes on each side (2 on the 
left and 3 on the right as depicted in Fig. lb, for example). These may be localized 
atomic orbitals or band modes of a solid crystal. (Small letters / and r are used here 
just to distinguish the single from the double junction.) Putting the left and right side 
together, scattering is determined by some complex amplitudes Sij gathered into matrix 
S (Fig.lc). The case that some modes may not couple is included. Some entries Sij 
may be zero. Reverse scattering is given by the adjoint matrix (Fig. Id). Not every 
matrix with as many lines as there are modes on one side and as many columns as 
there are modes on the other side can be a scattering matrix S, though. Probability 
conservation sets upper limits on the entries. A mandatory limit is, of course, that the 
absolute value of each entry be lower than or equal to 1, \sij\ < 1 for all i,j. But then 
one could think at first glance that it would be sufficient to demand that a particle 
occupying a pure mode on one side after being scattered to the other side should not 
cause a probability exceeding 1 there, that is < 1 and then the other way 

round J2j kjjP ^ 1 every j. An example showing that this is insufficient will be 
given in section 4. For independent transport channels we want an input on the left to 
come back into the same eigenmode it came from there after having been scattered to 
the right and back. All further multiple reflections will then stay in this mode. The 
channels will therefore be determined by diagonalizing S'^S. By construction S'^S is a 
hermitian matrix and thus has real eigenvalues. To ensure that the total probability 
does not increase it is required that any superposition of modes on the left described by 
a normalized distribution vector when scattered to the right and back is projected onto 
a mode distribution with total probability less than or at most 1 again on the left. The 
same, of course, has to hold true for starting with a normalized vector on the right and 
regarding the probability returning from one backreflection to the left. We shall show 
that the prohibition of probability creation in forth- and backreflections is equivalent to 
all eigenvalues of S'^S being less or equal to one. "All eigenvalues of S'^S < 1" is not 
a property to deduce of a scattering matrix S, but the definition to put for calling a 
matrix 5* a scattering matrix. 

In contrast to the usual scattering formalism [U [5] we do not distinguish incoming 
and outgoing modes here. Furthermore instead of two transmission and two reflection 
quadrants our scattering matrix S only consists of one transmission block and only 

I The model later will, however, assume that each mode has a constant density of states around 
the Fermi energy in a range a few times the equivalent of the voltage applied over the contact or the 
correspoding bulk superconductor quasiparticle density of states according, for example, to BCS theory. 
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Figure 1. (a) Point contact embedded in an electric circuit, (b) Modes left and right 
that can arbitrarily couple, (c) Scattering matrix from left to right, (d) Scattering 
matrix from right to left. 



describes transmission from one side of the contact to the other, and therefore does 
not have to be square. The term modes or original modes is used for states that 
conveniently describe the left and right side of the contact, which can be taken from 
the situation before a contact is established. We avoid the word "channel" for these, 
which is often used synonymous with incoming and outgoing modes [3l [5] because it is 
implied differently in the experiment-related description of atomic contacts [7] . We call 
eigenmodes or new modes those linear combinations of original modes that constitute 
eigenvectors of S'^S or SS"^ . Eigenmodes on the left and the right side of a contact are 
associated to each other one to one and their connections are called transport channels 
or simply channels. 

Let aij denote the entries of the matrix (j^t^) such that each column of (|^ts) is 
a normalized eigenvector of S'^S. Eigenvectors associated to eigenvalues zero have to 
be kept, and in case of degenerate eigenvalues choose a set of orthogonal eigenvectors 
associated to them, (ji^ts) is a square matrix with the dimension the number of modes 
on the left. The orthonormality of the eigenvectors is expressed through 

J2 a*jiajk = Sik. (1) 
j 

Any normalized input state on the left given by a distribution vector b in the basis of the 
original modes can be converted into the basis of eigenvectors, where we shall call this 
same vector c. (I^ts) is invertable. b = (([^15)^ or bm = Y.j o.mjCj. From Y.m Kn^m = 1 it 
easily follows that also X]j c*Cj = 1 and vice versa. The following calculation shows the 
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above claimed equivalance about S being a scattering matrix: 

^ J2{s^s)*^b*{s^s),,b,<i 

i,j,k,l,m 

E i^l^il)* ^mf^im Cm ^ ^ 

i,l,m 

^ E4A^c^<l (2) 

m 

J2k{S^ S)ikO'km — ^mdim expresscs that the mth column of (j^ts) is the eigenvector 
corresponding to eigenvalue A^. Of course, = A^, because eigenvalues are real, and 
(1) has been used. 

— * 

Some normalized vector b or equivalcntly any normalized vector c can arbitrarily 
be chosen, c could especially be a vector with any one component equal to 1 and all 
other components 0. Then obviously A^ < 1 ^ |Am| < 1 for every m individually. To 
initialize the above calculation from the last line, suppose that every |Am| is smaller or 
equal to 1. Then, for a general c, (2) simply means weighing each term in the sum 
J2m c*mCm = 1 by a factor A^ < 1, which must give a result < 1. 

Having obtained the scattering amplitudes from a microscopic physical model 
should ensure S being a scattering matrix and all eigenvalues of S^S less or equal 
to 1. However, and also because of the need to make up number examples here, an 
easier criterion than calculating all eigenvalues of S^S would be helpful. 

Ei(^'^kr<i (3) 

is a sufficient condition, although not a necessary one. 

As a further important aspect one may wonder whether defining transport channels 
as (eigen-)modes from the right being backscattered only into themselves would have 
made a difference from having required this property for (eigen-)modes on the left. This 
would mean looking for the eigenvalues and eigenvectors of SS^. As SS^ may have a 
different dimension from S^S, one of these matrices may have more eigenvalues than the 
other one. The conjecture, that one might get a greater number of transport channels 
as well as channels with different transmissions with an ansatz looking at backscattcring 
from one side instead of from the other, however, is wrong. We shall now demonstrate 
that, if SS'' is of greater dimension than S'^S, SS"^ has all eigenvalues that S'^S has and 
all its remaining additional eigenvalues are zero. Note again that columns of (l^^ts) 
the eigenvectors of S^S: 

{s^sr{ii^,,r = {\,sisri^s^sr (4) 

(A^ts) is the diagonal matrix containing the eigenvalues of S'^S in the order of the 
columns of (i^t^)- To get each column vector of {^tg) multiplied by the respective A^, 
(Agts) has to be multiplied to (l^^ts) from the right. S'^S, (l^fets) ^-^^ (-^sts) are all square 
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matrices with dimension the number of modes ni on the left side. Now multiply (4) by 
S from the left and as matrix multiplication is associative, we can view that as 

{ss^rs-^"^{ify,sr = ^"'^"'(ifet5)"'(A5t5)"' (5) 

(5) tells us that the columns of 'S'(|:^ts) are eigenvectors of SS^ - though not necessarily 
normalized - with the entries of (Xsts) ^ associated eigenvalues. SS^ is a square matrix 
with dimension the number of modes rir on the right side. 5" as well as 'S'(|:^t5) has 
rir lines and rii columns. Single upper indices in (4) and (5) denote the dimension of 
a square matrix, double upper indices the line and column number of a rectangular 
matrix. Every eigenvalue of S'^'S is an eigenvalue of SS\ too. However, we have 
not yet specified which number of modes, ni or n^, is the greater. For equal mode 
numbers ni = rir we have just proven that the sets of eigenvalues of S'^S and SS'^ are 
identical. Now firstly suppose that rii < rir- Then (5) is a statement about ni out of the 
rir eigenvalues and eigenvectors of SS\ and gives no information about the other rir-ni. 
Secondly suppose that rii > rir- As a preparation for the following argument we show 
that the scalar product of any two different columns of ^'(l^ts) vanishes. Elements of 
(4) multiplied by (frgts) i^om the left give exactly those scalar products, (ef-^ts) is the 
matrix with the complex conjugates of the eigenvectors of S^S as lines, or (Ij^ts)^- '^ij 
as before refers to the entries of (|^ts)- 

i,j,m 

^ml^Li^ij^jk = XI Knl^ka-mk = ^ikXk (6) 

i,j,m 

S'djgts) has rii columns, each of which is an n^-vector. In an n^-dimensional vector space 
there can, however, only be rir non-vanishing mutually orthogonal vectors. To fulfill (6) 
the remaining rii-rir column vectors must be the zero vector. With more modes on the 
left than the right there arc necessarily some modes or linear combinations of modes on 
the left that do not get transmitted through the junction at all. Now in (6) choose k = I 
and k one of those column numbers for which J2j Sijajk is zero for all i. Then it follows 
that the corresponding A^, which is an eigenvalue of S^S, is zero as well. For the single 
junction diagonalizing S'^S or SS'^ will lead to the same ensemble of transport channels. 
For illustration an example with numbers is given in section 4. Regarding the first and 
last expression in the equation for Z = A;, (6) also tells us that each A^ represents the sum 
of some absolute values squared. Therefore by its special contruction S^S not only is a 
hermitian matrix with real eigenvalues, but all eigenvalues are even greater or equal to 
zero. This is essential because we shall identify the square roots of the A^ as phaseless 
transmission amplitudes t. 

3. Double junction 

Properties of a series of two junctions will initially be given in terms of two scattering 
matrices, one between reservoirs L(left) and I(island), and one between reservoirs I and 
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Figure 2. (a) Scattering matrix and (b) adjoint for the double junction. 



R(right) (Fig.3a)|§| Multiplying each scattering matrix by its adjoint and diagonalizing 
would give transport channel ensembles for each junction separately and represent the 
setup for incoherent coupling [10] across the island. There are multiple reflections in 
each junction, but a charge transported to the island not to go back through that same 
junction relaxes there, that is changes the island's electrostatic potential, however, looses 
the information which (eigen-)mode it had come into. Diagonalizing matrices for each 
junction separately will not lead to the same linear combinations of island modes as 
eigenmodes for each side. Setting up our formalism to determine transport channels 
with coherent coupling between two the junctions is not obvious, should just any of 
the three reservoirs L,I,R be chosen to demand that eigenmodes there associated to 
channels exclusively return to themselves if transmitted away from that reservoir and 
backscattered into it. Nevertheless, as a scattering process starting on the island after 
two scattering steps offers no other possibility than to have brought the charge carrier 
back to the island, the above requirement for the island provides a promising ansatz. 

Therefore given scattering matrices from I to L and from I to R, with by definition 
the same number of I-modes, are just joint together into one scattering matrix S by 
taking all rows together. S and S"^ are formally given with exemplary numbers of modes 
three on the left, four on the island and two on the right in Fig. 2 (Fig. 3 are drawn 
for different numbers). Now S'^S is diagonalized, again denoting eigenvalues as Am and 
eigenvectors as (lii^ts)- The column vectors from (j^t^) build linear combinations of the 
original island modes. Eigenvectors of SS\ however, given by ^d^ts)! ^ow combine 
modes from both leads L and R. A a/A^ associated to an eigenvector out of (|%ts) or the 

§ Still the respective reverse scattering matrices are the complex conjugate transposed ones. With 
an extended island one could think of a charge being backscattered immediately at the island edge 
or crossing the island, being backscattered at the rear edge and thus resonator-like interference. For 
coherent coupling between junctions we are thinking of the island longitudinal dimension between the 
contacts comparable to the electron coherence length. The transverse dimension may be larger. In 
contrast to a quantum dot with discrete levels, we want to regard a bulk-like island. The set of original 
modes on the island should be the same in both scattering matrices, and therefore consist of eigenstates 
of this bounded space. We thus set the premises for still having a continuous and quasi-constant DOS 
as a function of energy, however, no extra phase factors to scattering amplitudes complicating the 
calculation for charges going across the island. 
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corresponding vector from ^(l:]^!^) gives us the transmission amplitude of a "combined 
channel" between the island and the entire outside world consisting of both leads. If 
there have been original modes at L and nji at R, every column of :^7|='S'(|35t5) is 
a normalized vector of length ul+uji, whose first components represent coefficients 
for L-modes and whose last ur components are coefficients for R-modes. However, like 
for the single junction case, modes at L and R are independent {JJ] Supplying a charge 
carrier in an original mode on the left, for example, makes the transmission eigenmodes 
it contributes to be populated with percentages the absolute values squared of their 
coefficients in the respective column from -^^di^t^). Therefore the transmission of an 
above mentioned combined channel splits up into the probability percentages left and 
right side modes contribute to the respective eigenmode. For assigning transmission 
amplitudes we take the roots of these parts to be multiphed by ^/X^. So for the 
transmission amplitudes ti and tn for the left and right side channels of a pair of 
channels coherently linked together on the island (=" combined channel") we obtain 



ti = \/A, 

and 

tR = \/X 



\ 



•n-L 



El^(^(fet5))-P (7) 



nL+riR 



\ E I^(^(fet5))-P (8) 

This construction demonstrates that a series of point contacts enclosing a bulk-like island 
even with coherent coupling across the latter is not equivalent simply to a number of 

II With the single junction we usuaUy do not regard the mapping of the original modes on a side onto 
the eigenmodes for transmission any more. For the supply of charge carriers the original modes are 
considered to be incoherently populated with probability one up the Fermi energy, which automatically 
makes the eigenmodes be populated each with probability one on each side for they are normalized 
linear combinations. 
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effective transport channels from left to right. Paths from left to right rest a channel 
through the left junction connected to a channel through the right junction and these two 
parts can have different transmissions. However, we cannot expect to describe coherently 
coupled junctions by the channel ensembles each one would exhibit as a single junction. 
Furthermore, by construction there is no coherent cross coupling on the island between 
our newly found channels. One channel from the left can only be coherently coupled 
to one channel from the right and vice versa (Fig.Sb). This finding elegantly solves 
the problem of generalizing a Green's functions method for modelling current-voltage 
characteristics in the coherent case [IT] to several channels per junction. The algorithm 
to calculate changing rates for the island charge is only needed for a single channel per 
junction. Contributions from all pairs of channels can then simply be added in classical 
rate equations. There may also be channels in a junction not coherently coupled to a 
partner channel in the other junction (or pairs with transmission amplitude zero in one 
half). Their island charging rates are even easier to calculate |10j. We should only expect 
such unpaired channels, however, if in the original pattern there are modes coupled 
across one junction that are totally decoupled from all others involved in a network of 
couplings over both junctions, like L3 and 15 in Fig.3a. 14 is not directly coupled to 
L, however indirectly via R3 and 13 as well as other paths, and will thus contribute to 
eigenmodes involving connections to both leads. Like in the single junction case, there 
may be specific linear combinations of original modes at every one of the three sites that 
are not transmitted at all. 

For a consistent description of the double junction with coherent coupling across the 
island, it has to be required that |S'^S'6| < 1 with the normalized vector b representing 
any linear combination of original modes on the island and 5* and 5*^ now the big matrices 
from Fig.2. The prohibition of probability creation here also implies that all eigenvalues 
Xm of S'^S have to be less than or at most one, even if the transmissions assigned to the 
channels in each junction further get multiplied with the parts by which a new combined 
island mode is transferred to each side. The new left and right modes at the ends of a 
pair of channels in fact constitute only parts of new combined-lead modes. The latter, 
however, form a unique set of eigenmodes. All these considerations about the double 
junction will be illustrated by number examples in the next section. 

We shall here present some further more marginal considerations on probability 
conservation or prohibiting probability creation. At the beginning of our analysis 
we regarded the system as described by an orthogonal set of modes on the left, an 
orthogonal set on the island and an orthogonal set on the right. Of course, with the 
contact established overlaps exist between modes in a lead and modes on the island. 
These determine the transition amplitudes put into the original scattering matrix. The 
overlaps should, however, only go so far that parts of an island wave function that 
couple to left wave functions rest well separated from parts that couple to the right. 
Then the sum of absolute squares of all scattering matrix elements out of an island 
mode will in no case exceed unity. The orthogonality of the modes in one lead ensures 
that the sum of absolute squares of all couplings to one certain island mode is limited 
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(c) 



by the norm one of that mode. Taking just any overlaps that could have been set up as 
single junction scattering amplitudes does not automatically guarantee that the norm 
of any column vector of the big matrix 5*, which contains scattering to the left and 
right modes, is less or equal to one. If the intersection of the overlaps of an island wave 
function with left wave functions and that with right wave functions is not empty, there 
would be a direct overlap between left and right wave functions, which will give rise to 
a direct transport channel from left to right. Such a situation could be realizable, but 
our model does not cover that. Nevertheless, suppose that two junctions in series have 
been investigated with incoherent coupling [10]. In the incoherent model changes of 
the island charge and thus potential are felt by both contacts, but otherwise the latter 
are independent and behave like a single junction each. With incoherent coupling the 
regarded set of island modes need not even be the same concerning transport through 
the left or the right junction, but here suppose it is. To compare with the incoherent 
behaviour, one might want to estimate the system's properties if the same junctions 
get linked together coherently (without allowing for direct left-to-right channels) and 
therefore just put together two single-junction scattering matrices Sli and Sri into one 
like in Fig. 2a. Where this should result in a column norm of S exceeding unity, in this 
situation it makes sense to renormalize scattering amplitudes in the way that of the 
island mode the part overlapping with both left and right wave functions is scattered to 
either side with probabilities proportional to the ratio of absolute squares of the single- 
junction scattering amplitudes. In any case put-together scattering matrices S should 
be checked for no column norm exceeding unity prior to solving for transport channels. 
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4. Examples with numbers 



(9) 



Let us illustrate that the simple criterion of norms of lines or columns being less or equal 
to 1 is not sufficient to ensure that a matrix is a scattering matrix. For simphcity an 
example with purely real numbers is given. For the single junction let be 

^ _ /0.3 0.5 0.5 0.6 \ 
"1^0.1 0.4 0.8 0.2 J 

supposing two modes on the left and four on the right. Then 

0.3^ + 0.5^ + 0.5^ + 0.6^ = 0.95, 0.1^ + 0.4^ + 0.8^ + 0.2^ = 0.85, 

0.41, 
0.40 



O.S^ + O.l^ 



0.10, 0.52 + 0.42 



0.52 + 0.8^ 



0.89, 0.6^ + 0.22 



(10) 





(11) 



(12) 



are all less than 1. 

/0.3 0.1\ 
0.5 0.4 
0.5 0.8 
V0.6 0.2/ 

The eigenvalues of this matrix, namely the solutions of 

(0.95 - A)(0.85 - A) - 0.75 ■ 0.75 = 0, 

however, are Ai=0.148 and A2=1.652 (rounded to 3 digits), of which A2 clearly is greater 
than 1. (3) is not fulfilled: 

0.952 ^ Q_752 ^ Q_752 ^ Q_g52 ^ 2.75 > 1 (13) 

Let us now regard 

/ 0.15 + 0i 0-0.50i \ 
S= 0.10-0.15i 0.10-0.05i , (14) 
V0.15-0.30i 0.05 + 0i / 
which is a scattering matrix, but otherwise its complex entries are randomly chosen. As 
S is 3x2, there are two original modes on the left and three on the right. 

S^S=( ^'^^^^ 0.025 - 0.05A 

V 0.025 + 0.05i 0.265 / ^ ^ 

and the eigenvalues of this 2x2 square matrix are Ai=0.1421 and A2=0.2904. 

/ 0.2725 0.04 - 0.0275i 0.0225 + 0.02i \ 

SS^ = 0.04 + 0.0275i 0.045 0.065 + 0.005i , (16) 

V0.0225 - 0.02i 0.065 - 0.005i 0.115 / 

is a 3x3 square matrix with eigenvalues Ai=0.1421, A2=0.2904 and A3=0 (a numerical 
calculation gives 2.7 • 10~^^). To check (3) for both S^S and SS^ we calculate 

0.1675^ + 2 



10.025 - 0.05i|2 + 0.2652 = 
0.27252 + 0.0452 + 0.1152 + 2 ■ |0.04 - 0.0275i|2 + 2 ■ |0.0225 + 0.02i| 
+ 2 • 10.065 + 0.005^|2 = 0.105 < 1 



:i7) 
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By the way, the same value is obtained as 
0.2904^ + 0.1421^ = 0.105 



(18) 



With complex numbers even normalized eigenvectors are only determined up to a phase 
common to all components. The numerical routine we used to evaluate them gave them 
with real last component. For S^S they are 

/ 0.407 - 0.814i\ /0.185 - 0.370A 

V -0.414 + 0i J \ 0.910 + 0i ) 

associated to Ai and A2, respectively. Eigenvectors of SS^ are 

-0.085 + 0.123A / 0.277 + 0.020i \ / 0.881 + 0.354i " 



(19) 



0.859 - OmSi 
-0.490 + Oi 



-0.457 + 0.043i 
-0.844 + Oi 



0.162 + 0.161i I (20) 
V 0.218 + 0i 
By multiplying (5) by 



associated to eigenvalues zero, Ai and A2 in that order. 
(ar^i^)S'^ from the left, one deduces that the factor needed to normalize a column 
vector of ^d^t^)) which we showed to be an eigenvector of SS\ is with A^ the 

corresponding eigenvalue. We shall now check explicitly the relation between the two 
sets of eigenvectors in our example, and especially verify that for unequal mode numbers 
left and right, the transformation in one direction produces some zero vector. 



V0.1421 



1 

V0.2904 



0.15 + 0i 
0.1 -0.15i 
0.15 - 0.3i 
/ 0.15 + 0i 
0.1-0.15i 
V0.15-0.3i 



- 0.5i 
0.1-0.05i 
0.05 + Oi 
- 0.5i 
0.1 - 0.05i 
0.05 + Oi 




0.407 - 0.814i^ 
-0.414 + 0i y 

0.185 - 0.372' 
0.91 + 0i 



^0.875i 



,-1.9i 



0.277 + 0.02?: 
-0.457 + 0.043* 
-0.844 + Oi 
/ 0.881 + 0.354i\ 

0.162 + 0.161i 
V 0.218 / 



(21) 



0.15 + 0i 
+ 0.5i 



0.1 + 0.15i 
0.1 + 0.005i 



0.15 + 0.3i 
0.05 + Oi 



0.15 + 0i 
+ 0.5i 



V0.1421 

0.1 + 0.15i 
0.1 + 0.005i 



0.15 + 0.3i 
0.05 + Oi 



0.15 -Oi 
+ 0.5i 



V0.2904 

0.1 + 0.15i 
0.1 + 0.005i 



0.15 + 0.3i 
0.05 - Oi 



0.277 + 0.02i 
-0.457 + 0.043i 
-0.844 + Oi 

0.881 + 0.354i~ 
0.162 + 0.161i 
0.218 + 0i 
-0.085 + 0.123i' 
0.859 - 0.003i 
-0.49 + Oi 



-0.875* 



0.407- 0.814i 
-0.414 + Oi 



A.9i 



0.185 - 0.37i 
0.91 + Oi 



1 ■ 10-^ + Oi 
5 • 10-5 + 1.5 • 10- 



(22) 



The last is the zero vector within the accuracy taking into account no more than 
the written digits. It cannot be normalized, of course, hence there is no factor 
in the last line. 
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Entries of scattering matrices are in general complex, but for simplicity we give 
examples with real numbers for the double junction, which is totally sufficient to show 
the important aspects of the calculation. To better keep track of which scattering 
amplitudes belong to the left and the right junction, S is noted in table form like in 
Fig. 2. Our first example has four modes on the island, two in each lead and no vanishing 
scattering amplitudes. 



S: 


h 


h 


h 


h 


u 


0.17 


0.28 


0.39 


0.06 


L2 


0.06 


0.22 


0.40 


0.11 


Ri 


0.28 


0.11 


0.33 


0.39 


R2 


0.06 


0.33 


0.11 


0.17 



(23) 



The sum of squares of all these 16 elements is 0.9921, which is less than 1 and should 
thus ensure that all eigenvalues will be less than or at most 1 . The matrix to diagonalize 
is 

/0.1145 0.1114 0.1893 0.1362\ 
0.1114 0.2478 0.2698 0.1400 

0.1893 0.2698 0.4331 0.2148 
V0.1362 0.1400 0.2148 0.1967/ 



S^S 



As eigenvalues and associated eigenvectors one finds 



(24) 



0.0080 



0.0518 



0.0959 



0.8364 



h 
h 
h 
h 



/ 0.8549 \ 
0.1327 
-0.2433 
V -0.4386/ 



/ 0.0728 \ 
-0.6389 
0.6469 
V -0.4100/ 



/ -0.3895 \ 
0.5854 
0.1872 
V -0.6859/ 



/ 0.3349 \ 
0.4812 
0.6981 

V 0.4111/ 



(25) 



S times the matrix made of these four column vectors returns a matrix consisting 
of the following four column vectors: 



Li 

Ri 
R2 



/ 0.0627 \ 
-0.0651 
0.0026 
V -0.0062/ 



/ 0.0612 \ 
0.0775 
0.0039 
V -0.2050/ 



/ 0.1296 \ 
0.1048 
-0.2504 
V 0.0738 / 



/ 0.4886 \ 
0.4504 
0.5374 

V 0.3256/ 



(26) 



(25) gives the normahzed eigenvectors of S^S given in the basis of the original 
island modes, which is hinted at by the labels li in front of their components. (26) lists 
eigenvectors of SS^ and their components refer to the basis of lead modes. Vectors in (26) 
are not yet normalized. To do so, divide by the square root of the above the respective 
column from (25). To fully link the number examples to the notation of the previous 
sections, remark that the vectors in (25) give the columns of {^tg) and those in (26) 
constitute ^'(l^ts)- Obviously the eigenchannel system consists of four pairs of channels, 
denoted as {tu, tRi), i = 1, . . . , 4. With the values from (26) we shall now work out the 
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left and right contributions in each. For example, the transmission amplitude of the left 
channel in the first pair is given by the part of the L-modes in the normalized eigenvector 
multiplied by the square root of the respective Am, which is the transmission amplitude 
of the island eigenmode to the entire outside world, tii = ^ 0-0627^^+ao65i^ . VO.0080. The 
reason to give non-normalized vectors in (26) was that explicit factors \/A^ drop out 
here. The following table hsts all left and right transmission amplitudes as determined 
by (7) and (8): 

i thi tm 



1 V0.06272 + 0.06512 ^ o.0904 

2 V0.06122 + 0.07752 ^ 0.0988 

3 V0.12962 + 0.10482 = 0.1667 

4 V0.48862 + 0.45042 = 0.6645 



V0.00262 + 0.00622 ^ o.0067 
VO.OO392 + 0.20502 = 0.2050 
VO.25042 + 0.07382 ^ 0.2610 
VO.53742 + 0.52562 = 0.6283 



(27) 



Remark that the sum of all ti squared and divided by the respective is 2 and 
the same holds for the t^. 



0.09042 0.09882 0.I6672 0.66452 

+ + + = 2.03 



0.0080 
O.OO672 



0.0518 
0.20502 



0.0959 
0.26102 



0.8364 
0.62832 



(28) 

_l_ |_ L — 1 9992 f29) 

0.0080 0.0518 0.0959 0.8364 ' ^ ^ 

Deviations from 2 here arc due to having limited accuracy to only four digits. (28) 
and (29) illustrate the reason and purpose of multiplying the \/A^ by the weights left 
and right sides have in the eigenvectors. Despite four channels ending in each lead, 
the maximum available supply of charges from each lead is reduced to the amount two 
original modes on each side could provide. We prefer to include these weight factors 
in tL and tu. Alternatively, the y/X^ could be called transmission amplitudes of the 
combined channels. However, then a channel would not connect to equal densities of 
states on the left and on the right, and weight factors for available states would have to 
enter effective transfer amplitudes, anyway. 

For comparison of the resulting channel transmissions let us now treat the system as 
incoherent, that is as two single junctions, with the same original scattering amplitudes. 



Sf. 


ii 


12 


13 


14 


h 


0.17 


0.28 


0.39 


0.06 


h 


0.06 


0.22 


0.40 


0.11 



SiS} = 



0.2630 
0.2344 



0.2344 
0.2241 



(30) 



SiSj has eigenvalues 0.0083 and 0.4786, which are also the non-vanishing eigenvalues 
of SjSi- They correspond to transmission amplitudes tn — ^0.0083 = 0.0911 and 
ti2 = V0.4786 = 0.6918. 
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ii 


12 


1.3 


14 


n 


0.28 


0.11 


0.33 


0.39 


r2 


0.06 


0.33 


0.11 


0.17 



SrSl 



0.3515 0.1557 
0.1557 0.1535 



(31) 



Eigenvalues of SrS} as well as of S^Sr are 0.0700 and 0.4370 and these translate 
into transmission amplitudes tri = \/0.0700 = 0.2646 and tr2 ~ -\/0.4370 = 0.6610. Like 
for the single junction the number of channels in each junction is equal to the smaller 
number of modes on one side of a junction. There are two channels per junction here. 
These may, however, have four different linear combinations of original island modes as 
new modes at their ends on the island. The transmission amplitudes are comparable in 
size to some of those found in the coherent case, mostly the greater ones for a junction, 
but there is no easy way of predicting the channel ensembles for the coherent double 
junction from those for both junctions taken as singles. 

In order to verify that modes only coupled across one junction and disconnected 
from a network extending over both result in channels through only one junction, we have 
done an analogous calculation to the above one starting with the following scattering 
matrix: 



S: 


h 


h 


h 


h 




0.17 


0.28 


0.39 







0.06 


0.22 


0.10 





Ri 


0.28 


0.11 


0.33 





R2 











0.17 



I 

1 

2 
3 
4 



0.0606 


0.1087 
0.6752 



0.0161 
0.17 
0.1497 
0.4204 



(32) 



This results in three eigenvectors the fourth component of which vanishes 
(eigenvalues 0.0039, 0.0342 and 0.6326) and one eigenvector of which the first three 
are zero (eigenvalue 0.0289). The resulting transmission amplitudes of channel pairs are 
listed in (32). The second line with tL2 = obviously represents a single channel only 
bridging the right junction. The input weight of left modes is 2, 

- (33) 



0.0606^ ^ 0.1087^ 
+ + 



0.0039 0.0342 0.6326 

The input weight of all right modes is also 2, 



2.01 . 



0.0161^ 0.17^ 
+ 



0.1497^ 0.4204^ ^ 

+ -TTTZTTTTT + = 2.001 



0.0039 ' 0.0289 ' 0.0342 ' 0.6326 ' ^^^^ 

where 1 falls to the single channel with iij2 = 0.17 and the others amount to 1 altogether. 

The following example is made up for a qualitative estimation on how a junction 
that exhibits a channel with medium transmission taken as stand-alone device and a 
tunnel junction that has many low-transmission channels would behave when combined 
with coherent coupling across the island between them. To make the system as coherent 
as possible except for paths through the tunnel junction, we assume that many modes 
from the right lead couple weakly to the same mode on the island. With one mode on 
the left, one on the island and modes i?i to i?io on the right let the original scattering 
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amplitudes be given by 





L 




/?2 




Rio 


I 


0.5 


0.01 


0.01 




0.01 



Then S^S 



0.251 and there is only one non-zero 



(35) 
v/(r25T = 



0.52 + 10 ■ 0.01^ 

0.501. The associated eigenvector of S'^S is (1), that is the island mode is the eigenmode. 
It translates into the combination 0.5L-|-0.01it!i-|-0.01ii!2 + - ■ .+0.01i?io in the lead basis. 
The system exhibits one channel pair with tL = VOI^ = 0.5 and tn = VlO • 0.01^ = 
0.032. Prom the right only the eigenmode (O.Oli^i + 0.01i?2 + . . . + 0.01i?io) /VlO • O.OP 
is transmitted, all other linear combinations of right modes are not. The tunnel-junction 
appears as one effective channel the transmission amplitude of which is enhanced over a 
single original path by the square root of the number of paths having been assumed with 
equal throughput. Another example with unequal numbers of original modes left and 
right will be contained in the next section. There is a tendency of left and right channel 
transmissions being paired ordered by size. This is no principal necessity, however. For 
example, the channel left with the largest (smallest) ti need not be coherently connected 
to the channel right with the largest (smallest) tR. 

The last example in this section is intended to illustrate how subtle the differences 
between a valid scattering matrix and a matrix violating probability conservation can be, 
especially in comparing double to single junctions. Condition (3) is a sufficient criterion, 
however, set up overcautiously, and thus not a necessary one. We have already seen that 
if we were to take the matrix from (9) for a single junction between L and I, for example. 



Si: 


ii 


12 


13 


14 


h: 


0.3 


0.5 


0.5 


0.6 


h: 


0.1 


0.4 


0.8 


0.2 



(36) 



although the norm of each line and each column vector is less than 1, this would 
result in an eigenvalue and thus a transmission amplitude for one of the two channels 
^2=V^=l-285 greater than one. Supposing that the first line gives scattering of four 
island modes to a single mode on the left, again for only a single junction between L 
and I, 

.'ir ii ir, io 

SiS}^ (0.95) (37) 





ii 


12 


13 


14 


I: 


0.3 


0.5 


0.5 


0.6 



then it is easily seen from SiSj that the only non- vanishing eigenvalue is 
A=0.95 corresponding to t=VO. 95=0. 975. Analogously supposing that the second line 
respresents scattering through a single junction, named as between island and right. 





ii 


12 


13 


14 


r: 


0.1 


0.4 


0.8 


0.2 



SrSl 



0.85 



(38) 



the eigenvalue and the transmission amplitude for the only channel here would be 
A = 0.85 and f=0.922. If now we were to take the matrix (9) as set up for the double 
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junction with one mode left and one right 



S: 


h 


h 


h 


h 


L: 


0.3 


0.5 


0.5 


0.6 


R: 


0.1 


0.4 


0.8 


0.2 



A: 

L : 
R : 



0.148 

0.6832 \ 
-0.7302 j 



1.652 

-0.7302 \ 
-0.6832 J 



(39) 



then from the corresponding eigenvectors of 5*5^ associated to Ai=0.148 and 
A2=1.652 we get as transmission amphtudes 



1 0.6832 ■ VoliS = 0.2628 

2 0.7302 ■ yL652 = 0.9385 



0.7302 • yol48 = 0.2809 
0.6832 ■ v/L652 = 0.8781 



1.0 
1.0 



(40) 
(41) 



and total mode weights left and right 

0.2628V0.148 + 0.9385Vl-652 
0.2809V0.148 + 0.878lVl-652 

Even though a A exceeds one, no t does. (36) is no allowed scattering matrix for a single 
junction, but (39) from the result for all t looks possible for a double junction. (We 
should nevertheless have rejected this S, because (13) tells us that one backscattering 
off the island of the two lead modes amounts to a weight of more than two, which 
indicates an overlap situation as in Fig. 4b.) These considerations show once more 
that the double junction has to be treated as a principally different situation and its 
quantitative properties cannot be deduced in a simple manner from those of the single 
junctions. 

For the algorithms presented here, that is for diagonalizing a matrix S'^S and 
determining its eigenvalues and even for splitting up such an eigenvalue into amplitudes 
ti and tn for a pair of channels in the double junction, it is of no importance whether the 
eigenvalues A are less or equal to one or not. You could be tempted to interpret a t > 1 
as summing configurationally degenerate modes in a reservoir that quantitatively have 
identical overlaps with all modes on the other side of a junction and together are well 
supplied with charge carriers with a weight two, three, etc. instead of one. Such modes 
should however be entered into 5* as several, albeit identical, lines or columns to ensure 
for the orthogonality of all modes and the premises of each being fed with weight one. 
Transmission amplitudes t greater than one are unphysical. Green's functions P, [THl [12] 
that represent transmissions renormalized for multiple reflections require single-hopping 
amplitudes less or equal to one per transport channel. 

Of course, being given scattering matrices for two individual junctions - taking 
reflection parts with - there is a straight forward rule how to combine two such matrices 

% The easiest case is T = </(! + t'^) for a channel in a single junction in the normal conducting case, 
which stems from a geometric series |13) . 
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into one describing the two junctions put in series (with what we call coherent coupling) 
[5J. So there has to be a reason why we reopen the problem from the point of view of 
eigenchannels and newly invent a way for describing the double junction. The fact that 
for a single junction from experimental results eigenchannel ensembles can be deduced 
whereas full scattering matrices cannot, does not lead to an irrefutable justification. 
Merely being given current- voltage characteristics, for the double junction there may 
even be some ambiguity in deciding whether transport is coherent or incoherent 
Like scattering matrices our concept with channel pairs primarily provides a theoretical 
ansatz for modelling. Convoluting scattering matrices of two junctions in series into 
one effective matrix once and for all [5j would be appropriate if the island between 
the junctions did not change its state through charge transport onto and off it (such a 
resulting matrix would mimic an effective single junction). However, in the situation 
we are interested in, prone to Coulomb blockade, the island potential changes with each 
charge transfer and thus the energy range of available electron or hole states does, too. 
Thus an all-inclusive scattering matrix or a renormalized transfer Green's function has 
to depend on all possible island charge states. Convoluting here scattering matrices 
or Green's functions of two junctions into effective global matrices or functions, with 
the further necessity to transform to Fourier space in order to account for energy- 
dependent densities of states, is practically impossible for cross-linked transfer paths as 
in Fig.3a, whereas this is feasible for paired channels as in Fig. 3b [12]. Besides that, 
the original transfer Green's functions approach to transport properties of especially a 
superconducting quantum point contact [9j was made without separating modes into 
ingoing and outcoming parts and with restricting scattering matrices to transmission 
parts (called hopping elements). An extension of that formalism to double junctions 
should keep these premises. 

5. Partial coherence 

With today's microstructering techniques, particularly with metals, it is possible to 
fabricate islands that are still large enough to exhibit a bulk-like continuous density of 
states, however already small enough to be sensitive to single-electron charging |14j . 
One may further presume a mixture of coherent and incoherent interaction between 
processes in both junctions. It shall be shown that such partially coherent transport 
can be described by finding a system of channels following exactly the same method we 
just presented for the coherent case. In fact, the system will look like a coherent one 
with channels of just little lower transmission than one set up with coherent coupling 
across the island only. 

Also the fully coherent case allows sequential (incoherent) transport. A coherent 
multiscattering process can end on the island letting an electron or hole having come 
in into any mode there relax into the reservoir of charge carriers equally populating all 

+ Island dimensions, for example for aluminum, arc typically from below one up to two micrometers 
across and some ten nanometers in height. 
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modes. Another process, that decharges the island and brings it back to its original 
potential, may then begin out of any mode on the island. 

To change from full coherence to partial coherence, we want to reduce coherence in 
transport across the island for a charge carrier coming into a mode there from a lead, 
but keep the coupling of that island mode to each lead the same, that is not decrease the 
transmissions of the junctions, an idea reminiscent of [15]. It can here be done by using 
the following picture. Replace the original island mode which was coupled to modes at 
L and R by a threefold of new modes: One that is coupled to L- and R-modes as before, 
but with all scattering amplitudes reduced by a factor of ^/p, one that is coupled only to 
the L-modes and one that is coupled only to the R-modes. As scattering amplitudes of 
the latter two new I-modes take the original amplitudes reduced by a factor of 
p is number between and 1 and can be interpreted as the probability that a charge 
carrier entering into the chosen island mode gets coherently transported across to the 
other junction. The sum of absolute squares of scattering out of the island mode stays 
the same as before. For L- or R-modes only their properties localized at a junction 
count. As compared to the single junction, in contrast to the island in the leads there 
is no newly introduced specific distance, over which it is important to know to what 
degree coherence is maintained or lost. For lead modes a splitting as done for island 
modes does not make sense. Of course, lead and island modes only coupled across one 
junction and in no way connected to a network coherently linking them to the other 
lead may always also be present from the beginning. 

For comparison with partial coherence we first elaborate the following 
fully coherent case with two modes on the island, two left and one right. 
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S: 


h 


h 


Li 


0.17 


0.30 


L2 


0.06 


0.40 


R 


0.28 


0.35 



0.1109 
0.1730 



0.1730 \ 
0.3725 J 



(42) 



S^S has eigenvalues with associated eigenvectors: 
0.0248 0.4586 



/ 0.8953 \ / 0.4455 \ 
1^-0.4455; 1^0.8953; 

Multiplying by S translates them into eigenvectors ol SS^: 



h 
h 



/0.17 
0.06 
VO.28 

/0.17 
0.06 
VO.28 



0.30^ 

0.40 

0.35, 

0.30^ 

0.40 

0.35. 



0.8953 \ 
-0.4455 ) 

■ 0.4455 \ _ 
0.8953 J ~ 



/ 0.0186 \ 

-0.1245 
V 0.0948 / 

' 0.3443 \ ^ 
0.3849 ^ 
.0.4381/ ^ 



R 



Li 
R 



(43) 



(44) 



(45) 



Prom left and right parts in these we read off the transmission amplitudes of two pairs 
of channels: 



1 \/0.01862 + 0.12452=0.1259 0.0948 

2 VO.34432 + 0.38452=0.5164 0.4381 

In this example remark that the input on the left having a chance to be transmitted 
amounts to less than 2 modes and that on the right to less than 1. 

0.125970.0248 + 0.516470.4586 = 1.22 (46) 

0.094870.0248 + 0.438l70.4586 = 0.78 (47) 

Together they have the weight of 2 modes. The island only offers 2 modes. Thus only 
a two-dimensional subspace of the three-dimensional vector space spanned by all lead 
modes can contribute to transmission through the system. 

Let us now suppose that the /2-mode is split into parts coherently transferred 
between the junctions and parts losing coherence on the island. The three modes 
replacing I2 are /2c for the coherent part and I2L and I2R for the parts coupling to 
the left and right lead only, respectively. We choose p — 0.5. In this special case, 
1 —p — 0.5, too. Compared to column I2 from S in (42) non-zero numbers in columns 
he, I2L and I2R here are multiplied by -\/(X5, for example, 0.30 • -\/(r5=0.212. 
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S: 


h 


he 


I2L 


I2R 


Li 


0.17 


0.212 


0.212 





L2 


0.06 


0.283 


0.283 





R 


0.28 


0.247 





0.247 



0.0692 \ 
0.0610 


0.0610/ 

(48) 

Because there is a lead mode less than island modes, one of the eigenvalues, noted 
below together with the corresponding eigenvectors, is zero. 



/0.1190 
0.1222 
0.0530 

V 0.0692 



0.1222 
0.1860 
0.1250 
0.0610 



0.0530 
0.1250 
0.1250 








0.0059 



0.1040 



0.3730 



h 

he 
hL 

hn 



/ \ 
-0.5774 
0.5774 
V 0.5774 / 



/ -0.7583 \ 
0.4085 
-0.0912 
V 0.4997 J 



/ -0.4353 \ 
0.0723 
0.6697 
V -0.5974/ 



/ 0.4852 \ 
0.7033 
0.4581 

V 0.2451/ 



(49) 



If now we multiply these by S we obtain the (non- normalized) eigenvectors of SS^ , 
noted again below the corresponding eigenvalues here: 







0.0059 



0.1040 



0.3730 



Li : /0\ /-0.0616\ / 0.0833 \ /0.3287\ (50) 
L2 : 0.0443 0.1839 0.3578 

R : \0j \ 0.0120 / V -0.2516/ V 0.3701/ 

The zero- vector here indicates that there is a combination of island modes h, he, 
hh and hn-i namely that given by the first vector in the above list of four-vectors, which 
is never transmitted. This state cannot exchange charge carriers with the lead modes. 
From the other three eigenvectors we read off the transmission amplitudes of three pairs 
of channels: 



1 V0.06162 + 0.04432=0.0759 0.0120 

2 VO.O8332 + 0.18392=0.2019 0.2516 

3 V0.32872 + 0.35782=0.4859 0.3701 

Here again, the maximum possible input left amounts to two modes, and that right 
to one: 

0.0759V0.0059 + 0.2019V0.1040 + 0.4859V0.3730 = 2 (51) 
0.012070.0059 + 0.251670.1040 + 0.370170.3730 = 1 (52) 

Channels are not at all the same as when letting the /2-mode being transported fully 
coherently across the island. The important finding is, however, that, even if by 
construction some of the original island modes only have a direct overlap with one lead, 
the eigenchannel ensemble does not contain channels merely bridging a single junction 
(if no paths are totally separated from the rest of the original network). The connection 
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of the coherently and incoherently transported island modes via the lead modes causes 
the system to appear as effectively consisting of fully coherently linked channel pairs. 
As already pointed out earlier, the system of channel pairs does not suppress sequential 
transport. The latter may well be found to provide the dominating contributions in 
current- volt age characteristics. The channels with the largest throughput left and right 
may form a pair, in principle, however, as well belong to different pairs 0- 

6. Conclusions 

It has been discussed how to determine the ensemble of transport channels for a series of 
two point contacts enclosing an island between them, which allows coherent transport 
across it. We find that transport channels are only coherently coupled together in 
pairs involving one channel per junction. This makes a Green's functions algorithm 
developed for single- channel junctions coherently coupled via an island [H] applicable to 
the general case of multichannel junctions. The channel ensembles describing coherently 
coupled contacts will differ from those describing the same contacts each taken as a 
stand-alone device. Partially coherent transport across the island can also be treated 
with the presented method of determining transport channels and will effectively look 
like coherent transport through more but less open channels. Circuits with more than 
two contacts in series or with short-cutting channels between leads, however, would be 
more complicated to handle. 
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